AD- A 136  432  ADAPTATION  OF  FLUX-CORRECTED  TRANSPORT  ALGORITHMS  FOR  1/| 
MODELING  OUSTY  FLOWS ( U )  NAVAL  RESEARCH  LAB  WASHINGTON 
DC  M  A  FRY  ET  AL .  20  DEC  83  NRL-MR-5230 
UNCLASSIFIED  F/G  19/4  NL 


/)£>-  A 136432 


UCuWlTY  CLASSIFICATION  OF  THIS  WAGE  '**•*  D*i* 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


I.  ACCIAiCnT'S  catalog  NUMACA 


4.  TITlC  (t>4  SoUttto) 

ADAPTATION  OF  FLUX-CORRECTED  TRANSPORT 
ALGORITHMS  FOR  MODELING  DUSTY  FLOWS 


t.  TVA*  04  ASAOAT  A  ACAIOO  COVCACO 

Interim  report  on  a  continuing 
NRL  problem. 


*•  ACAAOAMING  OAO.  ACAOAT  NUMACA 


1-  AUTMOAIw 


M.A.  Fry*  and  D.L.  Book 


ONTAACT  OA  GAAmT  NUMACAfa) 


N00014-82-C-2229 


»  ACAAOAMING  OAGANIZATION  KAMI  ANO  AOOACM 


Naval  Research  Laboratory 
Washington,  DC  20375 


CONTAOLLING  OAAiCC  n AMC  ANO  AOOASSI 


Defense  Nuclear  Agency 
Washington,  DC  20305 


10.  AAOGAAM  SLtMCNT.  AAOJCCT.  TASK 
AAtA  *  AOAK  UNIT  NUMACAS 


62715H;  44-0578-0-3 


12.  ACAOAT  OATS 

December  20, 1983 


•  I.  NUMACA  or  AAGCS 

23 


UONITOAIHG  aGCnCV  NAMC  a  AOOACSVJf  ««MM  Inm  Cmntmilnt  OWu)  IS.  SCCUAITV  CLASS.  tot  Oil*  wnj 


I*.  0IITAIAUTION  STATCMCnT  tot  0,1,  AaA art) 


UNCLASSIFIED 


Approved  for  public  release;  distribution  unlimited. 


Accession 


IT.  OlSTftltuTiON  ST ATCMCMT  (ol  th*  abtfraci  Mftrtd  In  Block  20,  It  Otttoront  from  Report) 


mmm 


1  •i-innounc 

■■  •  -  U'icn 


I*.  SUAALCMCN T AAV  NOTCS 

♦Present  address:  Science  Applications,  Inc.,  McLean,  VA  22102 
This  work  was  supported  by  the  Defense  Nuclear  Agency  under  Subtask  N99QAXAH, 
work  unit  00061,  work  unit  title  “Cloud  Calculations,”  and  under  Subtask  Y99QAXSG 
work  unit  00057.  work  unit  title  “Air blast 


19.  KCY  WO  AOS  (Confirm*  on  tovoro*  it  do  U  noeooomty  md  Identity  Or  ktock  ntmkor) 

Flux-corrected  transport  Airblast 

Finite-difference  algorithms  Height-of-burst 


20.WSITRACT  (Confirm*  m  toworoo  •  Ido  It  nocomomy  tmd  Identity  Of  ktock  iwrtwj 

Blast  wave  phenomena  include  reactive  and  two  phase  flows  resulting  from  the  motion  of 
chemical  explosion  products.  When  the  blast  wave  interacts  with  structural  surfaces  (external 
discontinuities),  multiple  reflections  and  refractions  occur  from  both  external  and  internal 
discontinuities.  The  most  recent  version  of  the  Flux-Corrected  Transport  (FCT)  convective- 
equation  solver  has  been  used  both  in  one  and  two  dimensions  to  simulate  chemical  explosive 
blast  waves  reflecting  from  planar  structures  for  yields  ranging  from  8  lbs  to  600  tons.  ^  ^ 

(Continues) 

DO  I  1473  SOITION  OA  I  NOV  «S  IS  OBSOLCTC 

s/n  oio  j-oi«-«soi 


IBCUAITV  CLASSIFICATION  OA  THIS  A  AOS  I 


'  i 
it 


t 


SCCUWITV  CUHl^lCATlQti  Q9  ThU  mgt  fWhmm  Of  toumn 


20  A0STHACT  fContiau+d) 

One  can  relate  the  strength  of  the  second  reflected  peak  to  the  sharpness  of  the  contact 
discontinuity,  and  thus  measure  the  capability  to  predict  all  the  salient  features  of  the 
blast  wave.  The  flow  patterns  obtained  reveal  four  different  vortices,  two  forward  and 
two  reversed.  Their  effect  on  the  motion  of  tracer  particles  has  been  studied  in  order  to 
determine  the  motion  of  (1)  HE  detonation  products  and  (2)  dust  scoured  up  from  the 
ground. _ 

'\ 


CONTENTS 


INTRODUCTION  . 

NUMERICAL  TREATMENT . 

600- TON  ANFO  EXPLOSION  AT  166  FT 

CONCLUSIONS . 

REFERENCES  . 

ACKNOWLEDGMENT  . 


ADAPTATION  OF  FLUX-CORRECTED  TRANSPORT  ALGORITHMS  FOR  MODELING  DUSTY  FLOWS 


M.  A.  Fry,  Science  Applications,  Inc. 
and 

D.  L.  Book,  Naval  Research  Laboratory 


INTRODUCTION 

In  this  paper  we  describe  a  series  of  calculations  carried  out  as  part  of  an 
ongoing  effort  aimed  at  studying  blast  wave  diffraction  effects  in  air.  The 
phenomena  of  chief  interest  to  us  include  velocity  fields,  particularly  those 
associated  with  the  toruses  (both  forward  and  reverse)  in  the  neighborhood  of  the 
rising  fireball,  and  the  distribution  of  dust  lifted  off  the  ground  by  the  winds 
and  the  structure  of  the  cloud  at  the  time  of  stabilization,  we  are  interested  in 
studying  the  nature  of  the  gas-dynamic  discontinuities  which  appear,  the  vortices 
(both  forward  and  reverse),  and  how  the  dust  content  of  the  ear  affects  the 
evolution  of  the  blast  wave. 

The  technique  we  have  employed  for  this  purpose  is  numerical  model  inq.  Oie- 
and  two-fluid  hydrocodes  based  on  the  Flux-Corrected  Transport  (PCT) 1  shock¬ 
capturing  techniques  have  been  used  to  simulate  airblast  phenomena  in  one  and  two 
dimensions.  PCT  refers  to  a  class  of  state-of-the-art  fluid  computational 
algorithms  developed  at  NRL  in  the  course  of  the  past  ten  years  with  supersonic 
gas-dynamic  applications  expressly  in  mind.  We  have  concentrated  on  modelina  the 
"Direct  Gburse"  event,  an  experiment  to  be  fielded  shortly  by  the  Defense  Agency: 
a  600-ton  ammonium  nitrate  +  fuel  oil  (ANPO)  charqe  is  detonated  at  a  height  of 
burst  (HOB)  of  166  ft.  The  results  are  most  conveniently  exhibited  in  terms  of 
velocity  vector  plots  and  tracer  particle  trajectories.  Examples  of  these  are 
presented  to  illustrate  our  results  and  conclusions. 

The  plan  of  the  paper  is  as  follows :  In  the  next  section  we  discuss  our 
numerical  techniques  and  validation  procedures.  In  Section  3  we  describe  the  600- 
ton  2D  HOB  calculations.  In  Section  4  we  summarize  our  conclusions  and  discuss 
their  domain  of  validity. 

Manuscript  approved  September  13,  1983. 
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NUMERICAL  TREATMENT 


As  described  by  Boris  and  Bode1 ,  FCT  is  a  finite-difference  technique  for 
solving  the  fluid  equations  in  problems  where  sharp  discontinuities  arise  (e.q. , 
shocks,  slip  surfaces  and  contact  surfaces).  It  modifies  the  linear  properties  of 
a  second-  (or  higher-)  order  algorithm  by  adding  a  diffusion  term  during  convec¬ 
tive  transport,  and  then  subtracting  it  out  "almost  everywhere"  in  the  antidiffu¬ 
sion  phase  of  each  time  step.  The  residual  diffusion  is  just  large  enough  to 
prevent  dispersive  ripples  from  arising  at  the  discontinuity,  thus  ensuring  that 
all  physically  positive  conserved  quantities  remain  positive.  FCT  captures  shocks 
accurately  over  a  wide  range  of  parameters.  No  information  about  the  number  or 
nature  of  the  surfaces  of  discontinuity  need  be  provided  prior  to  initiating  the 
calculation. 

The  FCT  routine  used  in  the  present  calculations,  called  JPBPCT  (an 
advanced  version  of  EEBPCT2),  consists  of  a  flexible,  general  transport  module 
which  solves  1-D  fluid  equations  in  Cartesian,  cylindrical,  or  spherical  geometry. 
It  provides  a  finite-difference  approximation  to  conservation  laws  in  the  general 
form: 


L:  f  4x3V  »  -  /  ♦(u-u  )*dA  +  f  r dA,  (1) 

c  6V(t)  SAftT-9  “  <5A(t) 

where  <j>  represents  the  mass,  momentum,  energy  or  species  mass  density  in  cell 
SV(t),  _u  and  Ug  represent  the  fluid  and  grid  velocities,  respectively,  and  t 
represents  the  pressure/work  terms.  This  formulation  allows  the  grid  to  slide 
with  respect  to  the  fluid  without  introducing  any  additional  numerical  diffusion. 
Thus,  knowing  where  the  features  of  greatest  interest  are  located,  one  can  con¬ 
centrate  fine  zones  where  they  will  resolve  these  features  most  effectively  as  the 
system  evolves. 

The  same  transport  routine  was  orployed  for  both  coordinate  directions  in  the 
2D  r-z  code  (called  FAST2D)  via  timestep  splitting.  A  Jones-Wilkins-Lee  (JWL) 
equation  of  state  (BOS)  was  used  for  the  detonation  products  and  a  real -air  BOS 
was  used  outside  the  HB-air  interface3 .  The  routine  was  written  in  the  form  of  a 
table  lookup,  using  interpolation  with  logarithms  to  the  base  16  computed  by  means 
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of  logical  shifts4.  By  thus  taking  account  of  the  architecture  of  the  machine  (in 
these  calculations,  a  32-bit-word  two-pipe  Texas  Instruments  ASC)  it  was  possible 
to  generate  very  efficient  vector  aode,  decreasing  the  time  required  for  SOS  cal¬ 
culations  to  a  snail  fraction  of  that  required  for  the  hydro.  The  EDS  specifies 
pressure  as  a  function  of  density  and  internal  energy.  In  mixed  cells  the 
combined  pressure  was  calculated  according  to  Dalton's  law. 

The  initial  conditions  were  taken  to  be  the  self-similar  flow  field  used  by 
Kuhl,  et  al.3,  corresponding  to  a  spherical  Qiapman-Jouquet  detonation  at  the  time 
the  detonation  wave  reaches  the  charge  radius  (Fig.  1).  This  was  propagated  with 
the  ID  spherical  code  until  the  detonation  front  attained  a  radius  just  atelier 
than  the  BOB,  at  which  time  the  solution  was  laid  down  on  the  2D  mesh  (Fig.  2). 

The  boundary  conditions  were  chosen  to  enforce  perfect  reflection  on  the 
ground  and  on  the  axis  of  symmetry  [(d^/^  k  *  where  <j>  *  0,0,^,  and  v^,  *  0] , 
where  "t"  and  "n"  denote  tangential  and  normal  components,  respectively,  and 
outflow  on  the  outer  and  the  top  boundaries  [(d^/dn)^  *  0,  where  <t>  *  p,p,vfc,vn]. 

For  the  2D  calculations  the  mesh  was  typically  200  x  100.  Fixed  qriddinq  was 
used  to  minimize  numerical  errors.  The  zone  sizes  were  2.1  m  x  2.1  m,  respec¬ 
tively.  For  the  late-time  calculations,  a  fixed  mesh  with  100  zones  in  the  radial 
and  200  zones  in  the  vertical  direction  was  used,  with  all  cells  of  dimension 
4.2  m  x  4.2  m. 

To  study  the  motion  of  dust  particles  in  the  flow  field  generated  by  the 
calculation,  the  simplest  model  describes  dusty  air  as  a  single  phase  with  density 
and  adiabatic  index  chosen  appropriately.  This  approach  ignores  the  properties 
associated  with  the  particulate  structure  of  the  dust  and  the  process  of  scouring 
by  which  the  dust  enters  the  air.  A  more  realistic  picture  results  if  we  treat 
the  dust  as  a  distinct  phase,  described  by  equations  of  mass,  momentum  and  energy 
conservation,  as  has  been  done  in  cne  dimension  by  Miura  and  Glass5.  The  dust 
equations  are  aoupled  to  those  describing  the  air  through  drag  and  heat  transfer 
terms. 
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It  has  been  pointed  out  by  Kohl6  that  in  such  a  treatment  a  dust  par¬ 
ticle  tends  to  become  entrained  in  the  prevailing  flow  over  a  distance  ~  10 3 
particle  diameters.  Thus  a  one-phase  description  is  satisfactory  whenever  par¬ 
ticle  sizes  are  at  least  a  thousand  times  smaller  than  the  anallest  length  scale 
in  the  hydrodynamics.  For  the  present  calculation,  this  scale  is  roughly  1  m,  so 
particles  smaller  than  1  rim  can  be  regarded  as  totally  entrained. 

When  the  mass  density  of  the  dust  component  is  small  compared  with  air 
(or  HE  product)  density,  a  further  simplification  is  possible.  Entrained  dust 
particles  can  be  followed  by  passive  advection.  That  is,  the  wind  fields  u*, 
Uy  are  taken  from  a  dust-free  hydrodynamics  calculation,  and  dust  is  advected  in 
these  fields  according  to 


dx 

3t 


(2) 


In  this  approximation  we  ignore  the  effect  of  the  momentum  and  energy  transfer  on 
the  air  phase.  The  same  approximation  can  be  used  for  larger  ( nonent rained ) 
particles  also,  provided  we  include  inertia  and  drag  by  using  the  force  law 

dV 

ra  3?  =  “  mge2  +  D(v  -  V)  ,  (3) 


where  V  is  the  particle  velocity,  g  is  the  acceleration  due  to  gravity,  and  D  is 
the  empirical  drag  coefficient  employed  by  Miura  and  Glass5. 

Equations  (2)  and  (3)  apply  best  in  the  limits  of  extremely  small  and 
extremely  large  particles,  respectively.  Although  they  restrict  the  scope  of  the 
calculation  (by  requiring  the  dust  content  to  be  anall),  they  have  the  computa¬ 
tional  advantage  of  allowing  us  to  obtain  time -dependent  dust  distributions  for 
many  different  choices  of  dust  size  spectrum  and  initial  distribution  from  a 
single  hydrodynamics  calculation. 
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600-TCN  ANPD  EXPLOSION  AT  166  FT 

The  yield  and  HOB  (600  tons  and  166  ft,  respectively)  in  the  calculation  were 
chosen  to  equal  the  values  used  in  the  Defense  Nuclear  Agency  Direct  Course 
experiment/  which  we  are  simulating.  The  Chapman -\Jouguet  parameters  used  to 
initialize  the  spherical  free-field  calculation  were  taken  to  be  those  for  the 
NH4ND3-fuel  oil  (ANPO)  mixture  used  as  the  explosive. 

Pigs.  3{a)-(c)  show  the  contours  of  HE  density  and  internal  energy  per  unit 
mass  and  the  velocity  arrow  plot  at  t*0,  just  before  the  reflection  at  ground  zero 
occurs.  Pigs.  3(d)-(f)  show  the  corresponding  plots  54  ms  later,  while  Figs. 
3(g)-(i)  show  than  after  245  ms.  Note  the  reflected  shock  proceeding  upward, 
reflecting  again  off  the  fireball,  and  propagating  back  in  a  downward  and  outward 
direction.  The  interaction  of  this  shock  with  the  radially  inward  flow  near  the 
ground  generates  the  reverse  vortex,  which  is  clearly  seen  in  Fig.  3(i).  Note 
also  the  positive  vortex  forming  near  the  top  of  the  grid  in  the  same  plot.  The 
latter  results  when  the  upward -propagat ing  reflected  shock  interacts  with  the 
radially  outward  flow  near  the  top  of  the  fireball;  it  is  not  produced  by  the 
buoyant  rise  of  the  fireball,  which  at  these  early  times  has  scarcely  begun. 

To  look  at  the  evolution  of  the  fireball  at  late  times,  we  reinitialized  on  a 
larger,  aoarser  grid,  representing  a  cylinder  400  m  in  radius  and  800  m  high.  The 
first  300  cycles  approximately  reproduce  the  early-time  results.  The  spherical 
shock  breaks  away  and  leaves  the  mesh.  The  flows  remaining  on  the  grid  are  now 
subsonic  everywhere.  Then  the  fireball  begins  to  rise  and  the  subsequent  develop¬ 
ment  is  due  to  the  combination  of  buoyant  rise  and  the  action  of  the  vortices  set 
up  by  the  early  shocks. 

Figure  4a,  which  shows  the  trajectories  of  passively  advected  tracer  par¬ 
ticles  over  the  time  interval  1.8  sec  to  3.97  sec,  displays  the  vortices  very 
clearly.  Fig.  4b  shows  the  particle  paths  for  the  time  interval  3.97  sec  to  7.34 
sec.  Notice  that  there  are  four  vortices  visible  in  the  plot:  two  positive  and 
two  reversed.  The  additional  small  vortices  are  apparently  a  consequence  of 
entrainment  by  the  major  ones.  As  far  as  we  know,  their  existence  has  not  been 
noted  previously. 
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It  is  clear  that  the  major  qualitative  features  persist  lonqest  in  the 
velocity  plots.  This  is  not  surprising,  as  the  circulation  patterns  represented 
by  the  vortices  have  essentially  infinite  lifetimes  in  the  absence  of  viscosity. 
We  have  run  out  to  stabilization  (not  shown  here)  and  have  found  that  these 
features  persist  in  the  velocity  plots  to  the  very  end.  At  all  times  t>0  the  oeak 
flow  velocity  in  the  HE  case  exceeds  that  in  the  comparable  point  source  solution. 
This  is  a  reflection  of  the  fact  that  the  Chapman -Jouguet  solution  at  a  radius  of 
10  m  has  a  pressure  peak  of  52  kbar,  vs  3  kbar  for  the  Sedov  solution  at  the  same 
radius.  The  means  that  the  former  starts  out  with  rruch  more  violent  notion,  i.e., 
fluid  velocities  an  order  of  magnitude  larger.  All  in  all,  in  many  respects  the 
HE  case  does  not  closely  resemble  a  point  source. 


CONCLUSIONS 

We  have  described  a  numerical  simulation  of  the  Direct  Cburse  Event.  The 
code,  gridding,  and  method  of  solution  are  the  same  in  the  two  calculations.  The 
following  conclusions  appear  to  be  among  the  most  important. 

(i)  The  flow  establishes  a  pattern  of  four  vortices,  two  forward  and  two 
reversed,  instead  of  one  of  each. 

(ii)  The  upper  vortex  forms  first,  followed  by  the  reverse  vortex  near  the 
axis  of  synmetry  and  the  ground.  Adjacent  HE  products  beqin  to  be 
entrained  into  a  positive  vortex  over  a  lcnqer  period  of  time,  several 
seconds. 

(iii)  Once  picked  up  (scoured)  off  the  ground,  dust  is  efficiently  trans¬ 
ported  upward  by  the  reverse  vortex  farther  from  the  axis. 

Phenomena  neglected  in  the  present  model  (e.g. ,  terrain,  conditions  in  the 
boundary  layer,  turbulence,  humidity,  etc.)  are  unlikely  to  alter  the  above  con¬ 
clusions,  which  mainly  depend  on  the  characteristics  of  the  solutions  in  the 
interior  of  the  mesh  and  over  long  periods  of  time.  Further  study  of  the  tracer 
particle  motions  we  have  calculated  is  likely  to  be  illuminating,  particularly 
when  we  begin  to  consider  the  evolution  of  various  initial  configurations  as  a 
function  of  particle  size.  In  closing,  it  is  appropriate  to  emphasize  the  far- 
reaching  significance  of  the  role  played  by  the  HE-air  interface  in  the  dynamics 
of  both  airblast  and  cloud  rise  phenomena,  and  the  importance  for  numerical 
simulation  of  correctly  tr  sting  this  interface. 
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Olcy  Attn  FCT 
Olcy  Attn  FCPR 


Chief 

Field  Command 
Defense  Nuclear  Agency 
Livermore  Division 
P  0  Box  808  1-317 
Livermore,  CA  94550 
Olcy  Attn  FCPRL 

Director 

Joint  Strat  TGT  Planning  Staff 
Offutt  AFB 
Omaha,  N3  68113 
Olcy  Attn  DOXT 
Oley  Attn  JLA 
Olcy  Attn  JLTW-2 
Olcy  Attn  NRI-STINFO  Library 
Olcy  Attn  EPFS 

Commandant 

Nato  School  (Shape) 

APO  New  York,  NY  09172 

Olcy  Attn  U.S.  Documents 
Officer 

Under  Secy  of  Def  for  Rsch  4 
Engrg 

Department  of  Defense 
Washington,  DC  20301 

Olcy  Attn  Strategic  6 

Space  Systems  (OS) 

Director 

BMD  Advanced  Technology  Center 
Department  of  the  Army 
P  0  Box  1500 
Huntsville,  AL  35807 
Olcy  Attn  1CRDABK-X 
Olcy  Attn  ATC-T 

Commander 

BMD  Systems  Command 
Department  of  the  Armv 
P  O  Box  1500 
Huntsville,  AL  35807 

Olcy  Attn  BMDSC-R  N  HURST 

Chief  of  Engineers 
Department  of  the  Army 
Forrestal  Building 
Washington,  DC  20314 
Oley  Attn  DAEN-RDM 
Olcy  Attn  DAES-MCE-D 
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Deputy  of  Chief  of  Staff  for 
0?S  i  Plans 
Departsent  of  the  Army 
Washington,  DC  20310 
01 cy  Aten  DAMO-NC 

Commander 

Barry  Diamond  Laboratories 
Department  of  the  Army 
2800  Powder  Mill  Road 
Adelphi,  MD  20783 

(CNWDI- Inner  Envelope: 

Attn:  DELHD-RBH) 

Oley  Attn  Chief  Div  20000 
Olcy  Attn  DELHD-I-TL  (Tech  Lib) 

Commander 

U.S.  Army  Armament  Material  Readiness 
Command 

Rock  Island,  IL  61202 
Olcy  Attn  MA  Library 

Director 

O.S.  Army  Ballistic  Research  Labs 
Aberdeen  Proving  Ground,  MD  21005 
Oley  Attn  DRDAR-BLE  J  Keefer 
02cy  Attn  DRDAR-TSB-S  (Tech  Lib) 
Olcy  Attn  DRDAR-BLT  V  Taylor 
Olcy  Attn  DRDAR-BLV 

Commander 

O.S,  Aray  Communications  Command 
Fort  Huachuca,  A Z  85613 

Olcy  Attn  Technical  Reference  Div 

Commander 

U.S.  Army  Concepts  Analysis  Agency 
8120  Woodmont  Avenue 
Eethesda,  MD  20014 

Olcy  Attn  M0CA-ADL  (Tech  Lib) 

Commander 

U.S.  Army  Engineer  Center 
Fort  Belvolr,  VA  22060 
Olcy  Attn  ATZA 

Division  Engineer 
U.S.  Army  Engineer  Div  Huntsville 
P  O  Box  1600,  Vest  Station 
Huntsville,  AL  35807 
Olcy  Attn  HNDED-SR 

Division  Engineer 

U.S.  Army  Engineer  Div  Ohio  River 

P  0  Box  1159 

Cincinnati,  OH  45201 

(Unclassified  Only) 

Oley  Attn  0RDAS-L  (Tech  Lib) 


Director 

U.S.  Army  Engr  Waterways  Exper 
Station 
P  0  Box  631 
Vicksburg,  MS  39180 

Olcy  Attn  WESSD  G  Jackson 
Olcy  Attn  WESSA  W  Flathau 
Olcy  Attn  J  Strange 
Olcy  Attn  WESSE  L  Ingram 
Olcy  Attn  Library 

Contender 

O.S.  Army  Foreign  Science  &  Tech  Ctr 
220  7th  Street,  ME 
Charlottesville,  VA  22901 
Olcy  Attn  DRXST-SD 

Commander 

U.S.  Army  Material  &  Mechanics 
Rsch  Ctr 

Watertown,  MA  02172 

(Address  CNWDI:  Attn: 

Document  Control  for) 

Olcy  Attn  Technical  Library 
Olcy  Attn  DRXMR  J  Mescall 
Olcy  Attn  DRXMR-TE  R  SHEA 

Commander 

U.S.  Army  Material  Dev  &  Readiness 
CMD 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Olcy  Attn  DRCDE-D  L  Flynn 
Olcy  Attn  DRXAM-TL  (Tech  Lib) 
Und  only 

Commander 

U.S.  Army  Missile  Command 
Redstone  Arsenal,  AL  35809 
Olcy  Attn  DRDMI-XS 
Olcy  Attn  RSIC 

Commander 

U.S.  Army  Mobility  Equip  R6D  CMD 
Fort  Belvolr,  MD  22060 

(CNWDI  to  Army  Mat  Dev 
&  Readiness  Command) 

Olcy  Attn  DRDME-WC  (Tech  Lib) 

Cotsaander 

U.S.  Army  Nuclear  6  Chemical  Agency 
7500  Baeklick  Road 
Building  2073 
Springfield,  VA  22150 

(desires  only  Icy  to  Library) 
Olcy  Attn  Library 


Deputy  of  Chief  of  Staff  for 
OPS  &  Plans 
Departaent  of  the  Array 
Washington,  DC  20310 
Olcy  Attn  DAMO-NC 

Coaaander 

Barry  Diamond  Laboratories 
Departaent  of  the  Array 
2800  Powder  Mill  Bead 
Adelphi,  MD  20783 

(CNWDI-Inner  Envelope: 

Attn:  DELHD-BSS) 

Olcy  Attn  Chief  Div  20000 
Olcy  Attn  DELHD-I-TL  (Tech  Lib) 

Coaaander 

U.S.  Amy  Armament  Material  Readiness 
Comand 

Bock  Island,  1L  61202 
Olcy  Attn  MA  Library 

Director 

U.S.  kray  Ballistic  Beseareh  Labs 
Aberdeen  Proving  Ground,  MD  21005 
Olcy  Attn  DBDAR-BLE  J  Keefer 

02cy  Attn  DBDAB.-TSB-S  (Tech  Lib) 

Olcy  Attn  DBDAR-BLT  V  Taylor 

Olcy  Attn  DRDAk-BLV 

Coaaander 

U.S.  Amy  Comunications  Coaaand 
Fort  Huachuca,  AZ  85613 

Olcy  Attn  Technical  Beferenee  Div 

Co  remand  er 

U.S.  Amy  Concepts  Analysis  Agency 
8120  Woodaont  Avenue 
Bethesda,  MD  20014 

Olcy  Attn  MOCA-ADL  (Tech  Lib) 

Commander 

U.S.  Array  Engineer  Center 
Fort  Belvoir,  VA  22060 
Olcy  Attn  ATZA 

Division  Engineer 
U.S.  Amy  Engineer  Div  Huntsville 
P  O  Box  1600,  Vest  Station 
Huntsville,  A L  35807 
Olcy  Attn  HNDED-SR 

Division  Engineer 

U.S.  Amy  Engineer  Div  Ohio  Blver 

P  0  Box  1159 

Cincinnati,  0M  45201 

(Unclassified  Only) 

Olcy  Attn  0BDAS-L  (Teeh  Lib) 


Director 

U.S.  Amy  Engr  Waterways  Exper 
Station 
P  0  Box  631 
Vicksburg,  MS  39180 

Olcy  Attn  WESSD  G  Jackson 
Olcy  Attn  WESSA  W  Flathau 
Olcy  Attn  J  Strange 
Olcy  Attn  WESSE  L  Ingres 
Olcy  Attn  Library 

Coaaander 

U.S.  Amy  Foreign  Science  &  Tech  Ctr 
220  7th  Street,  BE 
Charlottesville,  VA  22901 
Olcy  Attn  DRXST-SD 

Coaaander 

U.S.  Amy  Material  6  Mechanics 
Bach  Ctr 

Watertown,  MA  02172 

(Address  CNWDI:  Attn: 

Document  Control  for) 

Olcy  Attn  Technical  Library 
Olcy  Attn  DB3CMR  J  Meseall 
Olcy  Attn  DRXMR-TE  R  SHEA 

Coaaander 

U.S.  Amy  Material  Dev  &  Beediness 
CMD 

5001  Eisenhower  Avenue 
Alexandria,  VA  22333 

Olcy  Attn  DRCDE-D  L  Flynn 
Olcy  Attn  DBXAM-TL  (Tech  Lib) 
Und  only 

Coaaander 

U.S.  Amy  Missile  Comand 
Bedstone  Arsenal,  AL  35809 
Olcy  Attn  DRDMI-XS 
Olcy  Attn  RSIC 

Coaaander 

U.S.  Amy  Mobility  Equip  R&D  <20 
Fort  Belvoir,  MD  22060 

(CKWDX  to  Amy  Mat  Dev 
&  Beediness  Comend) 

Olcy  Attn  DBDME-WC  (Tech  Lib) 

Comender 

U.S.  Amy  Nuclear  &  Chemical  Agency 
7500  Baeklick  Bead 
Building  2073 
Springfield,  VA  22150 

(desires  only  ley  to  Library) 
Olcy  Attn  Library 
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Air  Tore*  Geophysics  Laboratory 
Hunscora  AFB,  MA  01731 

Olcy  ACCn  LWV  K  Thompson 

Air  Force  Institute  of  Technology 
Air  University 

Wright-Patterson  AFB,  OH  43433 

(Coes  not  desire  classified 
documents) 

Olcy  Attn  Library 

Headquarters 

Air  Force  Systems  Command 
Andrews  AFB 
Washington,  DC  20334 
Olcy  Attn  DLW 

Air  Force  Weapons  Laboratory,  AFSC 
Kir tland  AFB,  NM  87117 

Olcy  Attn  NTES-C  R  Henny 
Olcy  Attn  NTED  R  Matalucci 
Olcy  Attn  NTE  M  Plamondon 
Olcy  Attn  R  Guiee 
Olcy  Attn  SUL  W  Lee 
Qley  Attn  D£X 

Assistant  Chief  of  Staff 
Intelligence 

Department  of  the  Air  Force 
Washington,  DC  20330 
Olcy  Attn  IN 

Ballistic  Missle  Office/DE 
Air  Force  Systems  Command 
Norton  AFB,  CA  92409 

(Civil  Engineering) 

Olcy  Attn  DEB 

Ballistic  Missle  Office/MN 
Air  Force  Systems  Command 
Norton  AFB,  CA  92409 

(Minuteman)  MNNX 
Olcy  Attn  MNNXH  D  Gage 

Deputy  Chief  of  Staff 
Research,  Development,  &  ACC 
Department  of  the  Air  Force 
Washington,  DC  20330 
Olcy  Attn  AFRDQSM 

Deputy  Chief  of  Staff 
Logistics  &  Engineering 
Department  of  the  Air  Force 
Washington,  DC  20330 
Clcy  Attn  LEEE 


Commander 

Foreign  Technology  Division,  AFSC 
Wright-Patterson  AFB,  OH  45433 
Olcy  Attn  NIIS  Library 

Commander 

Rome  Air  Development  Center,  AFSC 
Criffiss  AFB,  NY  13441 
(Desires  no  CKWDX) 

Olcy  Attn  TSLD 

Stragetic  Air  Command/XPFS 
Department  of  the  Air  Force 
Offutt  AFB,  NB  68113 

Olcy  Attn  NRI-STINFO  Library 
Olcy  Attn  EPFS 

Department  of  Energy 
Albuquerque  Operations  Office 
?  0  Box  5400 
Albuquerque,  NM  87115 
Olcy  Attn  CTID 

Department  of  Energy 
Washington,  DC  20545 

Olcy  Attn  OMA/RD6T 

Department  of  Energy 
Nevada  Operations  Office 
P  0  Box  14100 
Us  Vegas,  NV  89114 

Olcy  Attn  Mail  &  Records 

for  Technical  Library 

Department  of  the  Interior 
Bureau  of  Mines 
Bldg.  20,  Denver  Federal  Ctr 
Denver,  CO  80225 

Olcy  Attn  Tech  Lib  (Unci  only) 

Director 

Federal  Emergency  Management  Agency 
1721  I  Street,  NW 
Washington,  DC  20472 

Olcy  Attn  Hazard  Eval  &  Vul 
Red  Div 

Aerospace  Corp. 

?  0  Box  92957 
Los  Angeles,  CA  90009 
Olcy  Attn  K  Mir els 
Olcy  Attn  Technical  Inf or 
Services 
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Agbabian  Associates 
250  N  Nash  Straac 
El  Segundo,  CA  90245 

Oley  Attn  M  Agbabiao 

Analytic  Servicea,  Zne. 

400  Army-Navy  Drive 
Arlington,  VA  22202 

Olcy  Attn  G  Hasaalbachar 

Applied  Theory,  Inc. 

1010  Westwood  Blvd. 

Los  Angelas,  CA  90024 

(2eys  if  unclass  or 
Icy  if  class) 

Olcy  Aten  J  Trulio 

Artec  Associates-,  Znc. 

26046  Eden  Landing  Bo ad 
Hayward,  CA  94545 
Oley  Attn  S  Gill 

Aveo  Research  &  Systems  Group 
201  Lowell  Street 
Wilmington,  HA  01887 

Olcy  Attn  Library  A830 

BDM  Corp. 

7915  Jones  Branch  Drlva 
McLean,  VA  22102 

Olcy  Attn  A  Levagnino 
Olcy  Attn  T  Neighbors 
Olcy  Attn  Corporate  Library 

BDM  Corp. 

P  0  Box  9274 
Albuquerque,  NM  87119 
Olcy  Attn  R  Hensley 

Boeing  Co. 

P  0  Box  3707 
Seattle,  WA  98124 

Olcy  Attn  M/S  42/37  R  Carlson 
Olcy  Attn  Aeroapeca  Library 

California  Research  4  Technology,  Zne. 
6269  Verlel  Avenue 
Woodland  Hills,  CA  91364 
Olcy  Actn  Library 
Olcy  Aten  K  Kreyenhagen 

California  Research  4  Teeh,  Zne. 

4049  First  Street 
Livermore,  CA  94550 
Olcy  Attn  D  Orphel 


Calspan  Corp. 

P  0  Box  400 
Buffalo,  NT  14225 
Olcy  Attn  Library 

Denver,  University  of 
Colorado  Seminary 
Denver  Research  Institute 
P  0  Box  10127 
Denver,  CO  80210 

(Only  ley  of  class  rpts) 
Oleu  Attn  Sec  Officer  for 
J  Wlsotskl 

EG4G  Washington  Analytical 
Services  Center,  Znc. 

P  0  Box  10218 
Albuquerque,  NM  87114 
Olcy  Attn  Librery 

Eric  H.  Wang 

Civil  Engineering  Rsch  Fee 
University  of  New  Mexico 
University  Station 
P  0  Box  25 

Albuquerque,  KM  87131 
Olcy  Attn  N  Baum 

Card,  Inc. 

7449  N  Natchez  Avenue 
Niles,  IL  60648 

Olcy  Attn  G  Neidhardt 
(Unci  only) 

General  Eleetrlc,  Co 
Space  Division 
Valley  Forge  Space  Center 
P  O  Box  8555 
Philadelphia,  PA  19101 
Oley  Attn  M  Bortner 

General  Electric  Co. -Tempo 
816  State  Street  (P  0  Drawer  QQ) 
Santa  Barbara,  CA  93102 
Olcy  Actn  DAS ZAC 

General  Resaarch  Corp. 

Santa  Barbara  Division 
P  0  Box  6770 
Santa  Barbara,  CA  93111 
Oley  Attn  B  Alexander 

Higgins,  Auld  Association 
2601  Wyoming  Blvd  NE 
Albuquerque,  NM  87112 
Oley  Attn  J  Bratton 
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Ill  Research  Institute 
10  w  35th  Street 
Chicago,  IL  60616 

Olcy  Attn  Documents  Library 
Olcy  Attn  R  Welch 
Olcy  Attn  M  Johnson 

Information  Science,  Xnc. 

123  W  Padre  Street 
Santa  Barbara,  CA  93105 
Olcy  Attn  U  Dudziak 

Institute  for  Defense  Analyses 
400  Army-Navy  Drive 
Arlington,  VA  22202 

Olcy  Attn  Classified  Library 

J  H  Wiggins  Co.,  Inc. 

1650  S  Pacific  Coast  Highway 
Redondo  Beach,  CA  90277 
Olcy  Attn  J  Collins 

Kaman  Avidyne 
83  Second  Street 
Northwest  Industrial  Park 
Burlington,  MA  01803 
Olcy  Attn  Library 
Olcy  Attn  E  Criscione 
Olcy  Attn  N  Hobbs 
Olcy  Attn  R  Ruetenik 

Raman  Sciences  Corp. 

P  0  Box  7463 

Colorado  Springs,  CO  80933 
Olcy  Attn  F  Shelton 
Olcy  Attn  Library 

Karaan  Sciences  Corp. 

Southern  California  Operations 
101  Continental  Blvd  Suite  855 
£1  Segundo,  CA  90245 
Olcy  Attn  D  Sachs 


Lockheed  Missiles  &  Space  Co.,  Inc. 

P  0  Box  504 
Sunnyvale,  CA  94086 

Olcy  Attn  TIC-Library 

Los  Alamos  National  Scientific  Lab. 

Mail  Station  5000 

P  0  Box  1663 

Los  Alamos,  NM  87545 

Olcy  Attn  MS  670/J  Hopkins 
Olcy  Attn  DOC  CON  for  M  Sanpford 
Olcy  Attn  DOC  CON  for  R  Whittaker 
Olcy  Attn  DOC  CON  for  MS  364 
(Class  Reports  Lib) 

Olcy  Attn  DOC  CON  for  C  Spillman 
Olcy  Attn  DOC  CON  for  A  Davis 
Olcy  Attn  DOC  CON  for  R  Bridwell 

Lovelace  Biomedical  6  Environmental 
Research  Institute,  Inc. 

P  0  Box  5890 
Albuquerque ,  NM  87115 

Olcy  Attn  R  Jones  (Uncles  only) 

Martin  Marietta  Corp. 

P  0  Box  5837 
Orlando,  FL  32855 

Olcy  Attn  C  Fotieo 

McDonnell  Douglas  Corp. 

5301  Boise  Avenue 
Huntington  Beach,  CA  92647 
Olcy  Attn  R  Halprin 

Merritt  Cases,  Inc. 

P  0  Box  1206 
Redlands,  CA  92373 

Olcy  Attn  J  Merritt 
Olcy  Attn  Library 


Lawrence  Livermore  National  Lab. 

P  0  Box  808 
Livermore,  CA  94550 

Olcy  Attn  DOC  CON  for  L-200  T  Butkovieh 

Olcy  Attn  DOC  CON  for  Tech  Inf or  Dept.  Library 

Olcy  Attn  DOC  CON  for  L-205  J  Hearst  (Class  L-203) 

Olcy  Attn  DOC  CON  for  L-90  D  Norris  (Class  L-504) 

Olcy  Attn  DOC  CON  for  L-437  R  Schock 

Olcy  Attn  DOC  CON  for  L-7  J  Kahn 

Olcy  Attn  DOC  CON  for  L-96  L  Woodruff  (Class  L-94) 

Olcy  Attn  DOC  CON  fro  L-90  R  Dong 


Meteorology  Research,  Inc. 

46-  V  Woodbury  Road 
Aitadena,  CA  91001 
Olcy  Attn  W  Craan 

Kachan  H.  Nevmark  Console 
Eng  Services 
3106A  Civil  Sag  Bldg. 

University  of  Illinois 
Urbana,  IL  61801 

Olcy  Aren  N  Nevsark 

Oak  Ridge  National  Lab. 

Nuclear  Division 
2-10  Lab  Records  Div 
P  0  Sox  X 

Oak  Ridge,  TN  37830 

Olcy  Attn  Civil  Def  Res  Proj 
Olcy  Aten  DOC  CON  for  Central 
Research  Library 

Pacifica  Technology 
?  O  Sox  148 
Dal  Mar,  CA  92014 
Olcy  Attn  C  Kent 
Olcy  Attn  R  3jork 

Physics  International  Co. 

2700  Merced  Street 
San  Leandro,  CA  94577 
Olcy  Attn  E  Moore 
Olcy  Attn  L  Sehraann 
Olcy  Aten  Technical  Library 
Olcy  Attn  F  Sauer 

R  &  D  Associates 
?  O  Sox  9695 

Marina  Del  Ray,  CA  9C291 

Olcy  Attn  Technical  Infer  Ctr 

Olcy  Attn  A  Latter 

Olcy  Attn.  A  Kuhl 

Olcy  Attn  J  Carpenter 

Olcy  Attn  C  MacDonald 

Olcy  Attn  R  Port 

Olcy  Attn  J  Levis 

Rand  Corp. 

1700  Main  Street 
Santa  Monica,  CA  90406 
Olcy  ATTN  C  Mow 

Sandia  Laboratories 
Livermore  Laboratory 
?  C  Sox  969 
Livermore,  CA  94550 

City  Attn  DOS  CON  for  Library 
6  Security  Class  Division 


Sandia  National  Laboratories 
P  0  Box  5800 
Albuquerque,  KM  87185 

(Accn  Mail  Services  Section 
for  Intended  Recipient) 

Olcy  Aten  Mail  Ser  Sec  W  Roherty- 

Olcy  Attn  Mail  Ser  Sec  3141 

Olcy  Attn  Mail  Ser  Sec  L  Vortsan 

Olcy  Attn  Mail  Ser  Sec  A  Chaban 

Olcy  Attn  Mail  Ser  Sec  L  Hill 

Science  Applications,  Inc. 

?  0  Box  2351 
La  Jolla,  CA  92038 

Olcy  Attn  Technical  Library 

Science  Applications,  Inc. 

1257  Tasttan  Drive 
Sunnyvale,  CA  94086 
Olcy  Attn  J  Dishon 

Science  Applications,  Inc. 

2450  Washington  Avenue 
San  Leabdro,  CA  94577 
Olcy  Attn  I  Maxwell 
Olcy  Attn  D  Bernstein 

Science  Applications,  Inc. 

P  0  Box  1303 
McLean,  VA  22102 

Olcy  Attn  M  Knasel 
Olcy  Attn  B  Chaobers  III 
Olcy  Attn  R  Sievers 
Olcy  Attn  J  Cockayne 

Southwest  Research  Institute 
P  0  Drawer  28510 
San  Antonio,  7X  78284 
Olcy  Attn  W  Baker 
Olcy  Attn  A  Wenzel 

SRI  International 
333  Ravensvood  Avenue 
Menlo  Park,  CA  94025 

Olcy  Accn  S  Abrahasaon 

Systems,  Science  &  Software,  Inc. 

P  0  Box  1620 
La  Jolla,  CA  92038 
Olcy  Attn  Library 
Olcy  Attn  D  Crine 
Olcy  Attn  T  P-iney 
Olcy  Attn  X  Pyatc 
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Teledyne  Srown  Engineering 
Custings  Research  Park 
Huntsville,  AL  35807 

Olcy  Attn  J  Ravenseraft 

Terra  Tek,  Inc. 

420  Wakara  Way 
Salt  Lake  City,  UT  84108 
Olcy  Attn  Library 
Olcy*  Attn  S  Green 
Olcy  Attn  A  Jones 

Tetra  Tech,  Inc. 

630  N  Rosenead  Blvd. 

Pasadena,  CA  91107 
Olcy  Attn  L  Hwang 
Olcy  Attn  Library 

TRW  Defense  &  Space  Sys  Group 
One  Space  Park 
Redondo  Beach,  CA  90278 
Olcy  Attn  I  Alber 
Olcy  Attn  Tech  Infor  Ctr 
02cy  Attn  N  Lipner 
Olcy  Attn  P  Bhutta 
Olcy  Attn  D  Baer 
Olcy  Attn  R  Plebuch 

TRW  Defense  &  Space  Sys  Group 
?  0  Sox  1310 

Sar.  Bernardino,  CA  92402 
Olcy  Attn  E  Wcng 
Olcy  Attn  P  Dai 

Universal  Analytics,  Inc. 

7740  V  Manchester  Blvd 
Playa  Del  Rey,  CA  90291 
Olcy  Attn  E  Field 

Veidlinger  Assoc.,  Consulting  Eng 
110  E  59th  Street 
New  York,  NY  10022 
Olcy  Attn  M  Baron 

Veidlinger  Assoc.,  Consulting  Eng 
3000  Sand  Kill  Road 
Menlo  Park,  CA  94025 

Olcy  Attn  J  Isenberg 

Westir.ghouse  Electric  Corp. 

Marine  Division 
Ker.dy  Avenue 
Sunnyvale,  CA  94088 
Olcy  Attn  V  Volz 

Code  4040  O.  Book  100  cys. 

Code  2628  20  cys. 
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